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Understanding thermal quenching of photoluminescence from first principles. 
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Understanding the physical mechanisms behind thermal effects in phosphors is crucial for white 
light-emitting diodes (WLEDs) applications, as thermal quenching of their photoluminescence might 
render them useless. The two chemically close Eu-doped Ba 3 SieOi 2 N 2 and BaaSieOgNi crystals are 
typical phosphors studied for WLEDs. The first one sustains efficient light emission at 100° C while 
the second one emits very little light at that temperature. Herein, we analyze from first principles 
their electronic structure and atomic geometry, before and after absorption/emission of light. Our 
results, in which the Eu-5d levels are obtained inside the band gap thanks to the removal of an 
electron from the 4f 7 shell, attributes the above-mentioned experimental difference to an auto¬ 
ionization model of the thermal quenching, based on the energy difference between Eusa and the 
conduction band minimum. For both Eu-doped phosphors, we identify the luminescent center, and 
we show that the atomic relaxation in their excited state is of crucial importance for a realistic 
description of the emission characteristics. 

PACS numbers: 65.40.-b,71.20.Eh, 71.20.-b, 71.70.Ej, 78.20.N- 


White light-emitting diodes (WLEDs) are seen as 
the most promising light source to replace incandescent 
and ultraviolet (UV) lamps. The phosphor-converted 
WLEDs rely on multiple layers of rare-earth (RE) doped 
photoluminescent materials called “phosphors” to down- 
convert (Stokes shift) the monochromatic UV or blue 
light from an electroluminescent diode into a broad- 
spectrum white light. The photoluminescent properties 
in RE-doped phosphors are most often based on radiative 
5d-4f transition of the dopant. Phenomenological models 
of the luminescence and Stokes shift of such phosphors 
are based on configurational diagrams mm ( see Fi g-0- 
In the case of RE-doped phosphors, one electron is ex¬ 
cited from the RE-4f band into the excited RE-5d band 
by a UV/blue photon coming from the diode (absorp¬ 
tion). Such excited RE-5d electron modifies the forces 
inside the material, changing the average atomic posi¬ 
tions from Vo to Xq. After dissipation of the energy as 
phonons, the excited electron spontaneously emits light 
and the crystal returns to its original ground-state geom¬ 
etry after yet another relaxation. 

Unfortunately for the WLEDs industry, other non- 
radiative mechanisms allow the RE-5d electron to relax 
to the ground state. For example, the electron can lose 
entirely its energy as heat by emitting phonons or transfer 
its energy to impurities or defects called “killer-center” 
via, for example, photo-ionization [3]. A phosphor can 
have high photoluminescence efficiency at low tempera¬ 
ture but little at working temperature. The existence of 
such thermal quenching seem to be sensitive to minute 
details of the host material and RE dopants. 

A striking example is given by two chemically 
close oxynitride materials: Ba 3 Si 60 i 2 N 2 :Eu and 

Ba 3 Si 6 0gN4:Eu (named generically BSON compounds). 


The two hosts materials have very similar electronic and 
structural properties 0] but when doped with Eu, they 
exhibit different luminescence properties and tempera¬ 
ture behavior. The green phosphor Ba3Si60i 2 N 2 :Eu E5K 
115] has an intensity that stays almost constant with 
temperature up to its working temperature, while 
BasSigOglXpEu is a bluish-green phosphor [HUMS! that 
exhibits an unfavorable decrease of the luminescence in¬ 
tensity with temperature. 

The experimental absorption and emission spectra as 
well as their temperature dependence have been mea¬ 
sured experimentally [5] and lead to a full width at half 
maximum (FWHM) of the two emission curves at 4 K 
of 1500 cm -1 and 1300 cm -1 for Ba3Si60i 2 N 2 :Eu and 
BasSigOglXpEu, respectively. These two FWHM are 
considered rather narrow and therefore it is expected 
that the luminescence in both compounds only comes 
from one luminescent center (one non-equivalent crystal¬ 
lographic position), despite the existence of several Ba 
inequivalent sites (two in Ba 3 Si60i 2 N 2 :Eu and three in 
BasSigOgN^Eu ), where substitution by Eu might occur. 

In this paper, we clarify the physical origin of the ther¬ 
mal quenching in these RE-doped oxynitrides, by com¬ 
plementing the available experimental data with first- 
principle results. Additionally, we confirm the unique¬ 
ness of non-equivalent luminescent sites that are present 
in both doped materials, and identify them. We believe 
our results are general and concern other RE-doped phos¬ 
phors. Due to the size of the cells needed to represent the 
doped materials, the Bethe-Salpeter equation (BSE) j2U] , 
presently the best technique to predict optical properties 
in solids, could not be straightforwardly used. We rely in¬ 
stead on a simpler mean-field approach, in which a fixed 
4f hole is introduced, either frozen in the core by means 
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Configurational diagrams 



FIG. 1: Comparison between the 4f-5d crossing and 

the Dorenbos auto-ionization model to explain the thermal 
quenching. 
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FIG. 2: Absorption and emission levels (eV) of the two 
BSON luminescent centers from experimental data or Go Wo 
calculations (in black) and present calculations (in red). The 
horizontal band positions at 0 K are shown in black, and 
at 300 K are shown in green (sometimes superimposed when 
there is no thermal effect). 


of a pseudopotential, or artificially forced inside the band 
gap so that the Euga electron is located inside the band 
gap. This approach allows us to draw, for the first time, a 
global and coherent picture of the absorption/emission, 
the geometry relaxation and the temperature effect on 
the luminescence. 

Models for the thermal quenching. The 5d-4f cross¬ 
ing decay model (Fig. 00 predicts that large atomic 
geometry changes lead to the crossing of the Eu, 5 d and 
Eu 4 / configurational energy curves. If the temperature 
is large enough, the excited electron of the 5d band over¬ 
come the activation energy barrier Ea and decays non- 
radiatively to the ground state. Based on this model, 
Blasse and Grabmaier [51] made predictions about the 
thermal quenching: stronger non-radiative decays occur 
when the zero phonon line energy difference E 2 — Eq is 
small and the parabola offset as well as the vibrational 
frequencies are important. As pointed out by Mikami 
and co-worker @1351123], the Ba 3 SigOi 2 N 2 :Eu compound 
has a smaller zero phonon line, equivalent parabola offset 
and larger phonon frequencies than BasSigOgN^Eu, in 
contradiction to the observed thermal quenching. 

In contrast, Dorenbos [25] proposed an auto-ionization 
model, see Fig. [T|j, arguing that the Eusa comes close 
to the conduction band minimum (CBM) from the host 
material well before the 5d-4f crossing point. We will 
show later that indeed this situation can be inferred from 
first-principles calculations. At working temperature, the 
Eu 5d electron is transferred to the conduction band, be¬ 
comes delocalized and mobile, and undergoes energy dis¬ 
sipation through other non-radiative mechanisms (trap¬ 
ping, killer centers such as defects, lattice vibration ...). 
Based on this model, the smaller the Eusa-CBM gap, the 
larger is the thermal quenching. Dorenbos’ empirical for¬ 
mula [22| CBM —Eu. 5 d = ^feV links the size of this gap 
to the quenching temperature To. 5 , at which the emission 
intensity has dropped by 50% from its low temperature 
value. From experimental thermal quenching data of the 
emission spectrum, Mikami et al. [25] could estimate the 


Eu 5d to CBM gap to be approximatively 0.6 and 0.2 eV 
for Ba 3 Si 60 i 2 N 2 :Eu and BagSigOglN^Eu, respectively. 

To validate the Dorenbos’ model, we have gathered 
experimental data and employed accurate many-body 
G 0 Wo ab-initio bandgaps of 6.88 eV and 6.45 eV for 
the two host materials This allowed us to build a 
global picture of the level positions for the absorption 
and emission in the two compounds, see Fig. [2] For the 
absorption, the 4f-5d gaps of 310 nm (4 eV) and 440 nm 
(2.8 eV) in Ba3Si60i 2 N 2 :Eu come from the peak and 
shoulder positions of the experimental absorption spec¬ 
trum [5]. The equivalent 4f-5d gap in BasSigOgN^Eu is 
located at 340 nm (3.65 eV). Similar consideration leads 
us to 5d-4f gaps of 535 nm (2.32 eV) and 480 nm (2.6 eV) 
for Ba 3 SigOi 2 N 2 :Eu and BagSigOgN^Eu, respectively. 
There is an additional blue shift of 0.06 eV at 300 K 
in the emission curve in the case of Ba 3 SigOi 2 N 2 :Eu [5]. 
Fig. [ 2 ] is then confirmed and completed (in red) by our 
first-principles results, as explained below. 

However, let us emphasize the meaning of the energy 
changes that are gathered in Fig. [2j At variance with 
standard density-functional theory (DFT) band struc¬ 
tures [23] or even GW quasi-particle band structures [5B], 
the levels inside the gap are not associated to the removal 
or addition of an electron but instead to neutral excita¬ 
tions. The Eu 4 f level is occupied by seven electrons in the 
ground state (half-filled shell). The absorption of a pho¬ 
ton leaves a hole in this Eu 4 f shell and the electron that 
is transferred in the 5cl lower levels or in the CBM feels 
the attractive potential of the hole, which is absent in 
traditional band structure methods. On the other hand, 
the band edges keep their usual significance: they weakly 
feel the localized 4f hole, since the corresponding electron 
(or hole) band edge states are delocalized. 

The ground state. The crystal structure of the two 
hosts, where the Ba atom has two inequivalent positions 
in Ba 3 SigOi 2 N 2 and three in Ba 3 SigOgN 4 , is presented 
in Ref. 3|. By replacing one Ba atom of the primitive 
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cell by an Eu one, the doping concentration is 33%, much 
higher than the 2-10% experimental doping [T5| ■ There¬ 
fore, we relied on doped supercells that are three times 
larger than the primitive ones, EuiBagSiigOggNg and 
EuiBagSiig027Ni2, leading to an Eu doping concentra¬ 
tion of 11%. More details on the five supercell models are 
given in the Supplemental Material at [URL will be in¬ 
serted by publisher]. The latter also contains a detailed 
description of the methodology and numerical parame¬ 
ters of the calculations, whose brief description is given 
hereafter. 

The ground state calculations using the generalized 
gradient approximation (GGA) (27, +U (U=7.5 eV and 
J=0.6 eV) [28] [29] lead to Eusd bands that are not local¬ 
ized inside the bandgap of the host material, in contra¬ 
diction to both the 5d-4f crossing decay model and the 
Dorenbos auto-ionization model. For BagSigO^l^Eu, 
this had already been noticed in Ref. m- 

Depletion of the ff shell. As already mentioned, the 
presence of Eu.sd states inside the band gap in Fig. [2] is 
due to the hole left by the Eu 4 f to Eu, 5 d transition, creat¬ 
ing an attractive center for the Eusd electron. The Eusd 
electron and the Eu 4 f hole are strongly correlated and 
their description would require a cumbersome treatment 
at the BSE level. To make the problem tractable, we use 
two simpler mean-field approaches that allow us to de¬ 
scribe the Coulomb effect of the hole while the Eusd state 
is occupied. In the first one, the pseudopotential core¬ 
hole approach (see for example Ref. 01]), a Af 6 shell is 
frozen in the pseudopotential, instead of the ground state 
4f 7 shell. In the second approach, the Eu 4 f hole is treated 
explicitly in the valence manifold through forced occu¬ 
pation numbers of the Eu 4 f orbitals (equal occupation 
numbers). This approach is called “constrained DFT” 
(CDFT) and has been applied successfully to a wide 
range of materials [32H35] , including rare-earth doped 
solid-state compounds, by Canning and coworkers, to 
study the luminescence in Ce-doped [HI 07] and Eu- 
doped [35] inorganic scintillators. In both approaches, 
the Coulomb interaction is correctly treated between the 
Eu core and the explicitly treated Eusd electron, as well 
as between the electrons and the other nuclei, allowing 
to include the most relevant physics in the computation 
of the forces, and thus yielding reasonable atomic ge¬ 
ometries. We have used the pseudopotential core-hole 
approach to relax the geometries from first principles. 

Self-interaction in the excited state. However, for the 
comparison of the 5d electronic level with the CBM, both 
approaches are quite deceiving, as they include a non¬ 
physical self-interaction of the excited electron in the 
Eusd and they poorly describe many-body effects in both 
Eusd and the CBM. The self-interaction overscreens the 
electronic properties and push the Eusd too high in en¬ 
ergy, comparatively to the CBM. 

In order to describe the electronic properties properly, 
a self-interaction free method has to be used. We found 


out that even the simple Kohn-Sham band structure ob¬ 
tained with the depleted Emu shell charge density, but 
without occupying the Eusd orbitals, nor the conduction 
band minimum has no (or little) self-interaction bias be¬ 
tween these states. 

Moreover, the use of GGA+U is not enough to local¬ 
ize the Eu 4 f states inside the bandgap when an electron 
is removed. An additional atomic potential energy shift 
of 0.3 Ha for Ba3Si60i2N 2 and 0.28 Ha for BaaSieOgNd 
has to be applied on the Eu 4 f states in order to reproduce 
the known experimental data shown in Fig. [2j Since this 
additional atomic potential is applied in the same way 
for the two sets, we can still discuss their relative differ¬ 
ence. We propose to name our hybrid method “depleted- 
shifted 4f approach” as we remove one electron and apply 
an additional atomic potential energy to those 4f states. 
An homogeneous negative background is also added to 
preserve charge neutrality. 

In order to validate it, we have performed tests, for 
the primitive cell only, based on the BSE and the GW 
methodology. In particular, we note that the BSE builds 
an effective Hamiltonian that takes into account two 
Feynman diagrams: the direct (attractive Coulomb inter¬ 
action) diagram and the exchange diagram [2D]. The at¬ 
tractive screened Coulomb interaction obtained from the 
BSE is accounted for to a large extend in the core-hole 
and depleted-shifted 4f techniques because the change 
of the electronic potential resulting from the hole in 
the pseudopotential is computed self-consistently within 
DFT (screened electron-hole interaction) and because all 
4f states are strongly localized. In contrast, the exchange 
term V xc within DFT is only a one-point function and 
therefore does not have the adequate form to reproduce 
the four-point kernel function of the BSE. To evaluate the 
crudeness of the mean-field approach, we have computed 
in a restricted space of valence and conduction states, the 
Bethe-Salpeter exchange term for the BaaSigO^l^Eu 
compounds. We have obtained that the effect of the ex¬ 
change could be as large as 0.25 eV in high-lying excitonic 
states, but the lowest transitions were nearly unaffected 
(0.02 eV). Also, we have performed GW calculations on 
the doped primitive cell, with the core-hole Eu pseudopo¬ 
tential. The relative changes, between the Eu 5 d and the 
conduction band minimum, compared to the Kohn-Sham 
eigenenergies, were much smaller than the GW correc¬ 
tions for the valence to conduction band gap of the host 
materials. The latter corrections were on the order of 
2.1-2.2 eV for the host materials [4j, while the GW cor¬ 
rection to the Eusd to CBM gap were only on the order 
of 0.2-0.3 eV. 

The 5d and conduction band levels in the case of 
neutral excitations. The 5d and conduction band lev¬ 
els, obtained with a mean-field 4f hole in the depleted- 
shifted 4f approach, were computed for the rhombohe- 
dral Eu-doped supercells with different geometries. The 
results with ground-state atomic positions (absorption) 
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FIG. 3: Spin-polarized (up in black and down in red) bandstructures corresponding to the absorption (ground-state relaxed 
atomic positions) and emission (excited-state relaxed positions), in a rhombohedral supercell. One of the nine Ba atoms is 
replaced by an Eu one at the I, II or III inequivalent Ba positions, leading to an effective Eu doping of 11%. The zero energy 
is positioned at the top of the valence band of the host materials and is not shown here. For the emission, the (f) and (h) cases 
are the luminescent ones whereas the 3 others are not due to the fact that their Eusa band is localized above the CBM of the 
host material. 


are shown in Figs. [3])a)-(e), while those associated with 
the emission process, computed at the relaxed atomic 
positions, are presented in Figs. [3](f)-(j). The differences 
between ground-state and excited geometries are large, 
and crucial to correctly describe the emission process: 
in all five cases, the absorption geometry leads to one 
(or two) Eusa level(s) below the CBM (dispersive band), 
while this only happens in two cases for the relaxed geom¬ 
etry. The contraction of the O cage around the Eu atom 
can be extremely important for some bonds, e.g. up to 
16% contraction of the Eu-0 bond length with respect 
to the undoped case. This level of additional contraction 
is well in line with the change of cation radii associated 
with the change of oxidation state for the Eu atom that 
we describe. Indeed, the experimental Eu cation radius 
goes from 131 pm for Eu 2+ to 108.7 pm for Eu 3+ [55] , 
From the analysis of the emission in Fig. [3] we con¬ 
clude that the two Eui substitutions in BagSiisO^Ng 
and BagSii 8 C> 27 Ni 2 are the only luminescent centers due 
to the position of their Eusd band below the CBM, as 
expected from the small FWHM. We found out that the 
two gaps are 0.52 eV and 0.37 eV for EuiBagSii8036N 6 
and EuiBasSii 8027 Ni 2 , respectively. 

For the two luminescent centers, the 4f states (not 
shown in figure [3]), went up by 2.5 eV and 1.5 eV, due 
to atomic relaxation. Although the absolute position of 
the 4f state depends on the value of the U and the 4f po¬ 
tential shift, the comparison made here is relevant since 
the only difference between the absorption and emission 
calculations are the atomic positions. To strengthen our 
analysis, we notice that the lowest conduction bands of 
the three non-luminescent centers are spin degenerate be¬ 
cause the delocalized states are far from the Eu atom and 
therefore do not react to its magnetic moment. We have 
also made an analysis of the change of 4f state energy 
with respect to the 5d-CBM gap, for intermediate ge¬ 


ometries. We observe that the 5d-CBM gap closes well 
before the 4f energy comes close to the 5d energy, sup¬ 
porting the Dorenbos model. 

The temperature dependence. Finally, we also studied 
the effect of temperature on the CBM of the two BSON 
host materials, using the static Allen-Heine-Cardona 
(AHC) theory prOl - Ho] . leading to a similar renormaliza¬ 
tion of -0.17 eV for both compounds at 300 K. At present, 
AHC calculations for the doped materials are out of reach 
computationally. Still, we do not expect large relative 
differences between the two materials, and in any case 
much smaller than the -0.17 eV renormalization of the 
CBM. 

Conclusion. All of the above mentioned results allow 
us to complete our schematic drawing of Fig. [2j yielding 
a global picture on which theory and experiment agree. 
From Fig. [3] we also identify the only active luminescent 
center in both Ba 3 Si 6 0i 2 N2 and Ba 3 Si 6 0gN 4 . The three 
other non-equivalent centers have their Eu, 5 d states com¬ 
pletely above the CBM level of the host material there¬ 
fore leading to non-radiative emissions. This confirms 
the experimental results that show narrow and well re¬ 
solved emission peaks for the two BSON compounds. 
Fig- @ has been obtained through the depleted-shifted 
4f scheme, whose validity was checked by comparing it 
to the Bethe-Salpeter theory in a simplified case. As 
seen in Figure [3j the relaxation of the Eu doped excited 
state is of crucial importance as there are major changes 
in the electronic bandstructure before and after Stokes 
shifts. We confirm the Dorenbos auto-ionization model: 
the Eusd-CBM gap is 0.09 eV larger in Ba 3 SigOi 2 N 2 :Eu 
than in Ba 3 Sie0 9 N4:Eu. This energy difference is enough 
to show significant thermal quenching behavior differ¬ 
ences Ul]. The techniques that we have used, for the 
ground-state, the excited state, and finally, for the study 
of the gap with the conduction band, should be widely 
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applicable for realistic models of phosphors at a much 
lower computational cost than the much more expensive 
BSE technique 
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